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"LXMlf I17a mNSn XN MIXES SUBSONIC AND 
SUPEESONIO FLOW OF OOMPEBSSIBLB FLUIDS 
By Hsue— shea Sslen 



It is well known that the vortlcity for any fluid element 
is constant if the fluid is non— viscous and the change of 
state of the fluid is isentropic. When a solid "body is placed 
in a uniform stream, the flow far ahead of the hody is irrota-i- 
tional. Then if the flow is further assumed to he isentropic, 
the Torticity will he zero over the whole field of flow, "In 
other words, the flow is irrotational, ■ For such flow over a 
solid "body, it is shown hy Theodorsen (reference 1) that tho 
solid hody experiences no resistance. If the fluid has a 
small viscosity, its effect will he limited in the "boundary 
layer over the solid hody and the T50dy will have a drag duo 
to the skin friction, tDhis type of essentially isentropic ir- 
rotational flow_ is genarally ohserved for a streamlined hody "" 
placed in a uniform streamy if the velocity of the stream is 
kept helov/' tho so-called 'fcritical speed," 

At the critical speed or rather at a certain value of the 
ratio of the velocity of the undisturhed flow and the corre- 
sponding velocity of sound, shock waves appear. This phenome- 
non is called the "compress ihility hurljlo, '* Along a shock 
wave, tho change of state of the fluid is no longer isentropic 
although still adiaha"bic, IJhis results in an increase in en- 
tropy of the fluid and generally introduces vorticity' In an '~ 
originally irrotational flow. The increase in entropy of the 
fluid is, of course, the consequence of changing part of tho 
mechanical energy into heat energy. In other words, tho part 
of fluid affected hy the shock wavo has a reduced mochanical 
energy. Therefore, with tho appearance of shock waves, tho 
wako of the streamline hody is very much widened, and tho 
drag increases drastically,. Purthcrmoro, the accompanying 
change in the pressure distrihution over the body changes the 
aerodynamic moment acting oji it and in the case of an airfoil '. 
decreases the lift force. 
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All these conseq.uences of the breakdown of isentropic 
irrotational flow are generally undesira'ble in applied aero- 
dynamics. Its occurrence should he delayed as much as possl— 
"ble hy modifying the shape or contour of the hody. However, 
such endeavor will he very much facilitated if the cause or 
the criterion for the breakdown can he found fi:;_st. 



OEIEEEION rOR THE BRIAEDOWN Off 
ISENTROPIC IRROTATIONAL FLOW 



Taylor and Sharman (reference 2) calculated the succes- 
sive approximations to the flow around an airfoil hy means of 
an electrolytic tank. They found that when the maximum veloc- 
ity in the flow reaches the local velocity of sound, tho con- 
vergence of the successive steps seems to break down. This 
fact led to the identification of critical speed or critical 
Maoh number with the, Mach number of the undisturbed flow for 
which the local velocity at some point roaches the local vel- 
ocity of sound. However, there is no mathematical proof for 
tho coincidence of the critical Mach number so defined and 
the breakdown of isentropic irrotational flow, furthormore, 
such a definition for critical Mach number implies that a 
transition from a velocity loss than that of sound, or sub- 
sonic velocity, to a velocity greater than that of sound, or 
supersonic velocity, does not occur in isentropic irrotational 
flow. On the other hand, Taylor (reference 3) and others 
found solutions for which such a transition ocetlF^, Turther— 
more, Binnie and Hooker (reference 4) have shown that at 
least for the case of spiral flow the method of successive 
approximation is a convergent one even for supersonic vel- 
ocities. With these facts in mind, i-t may be concluded that 
tho identification of critical speed with local supersonic 
velocity cannot bo correct, 

Taylor«s investigation on tho spiral flow (reference 3) 
indicates that there is a line in the flow field where the 
maximum velocity is reached and beyond which the flow cannot 
continue, Tollmien in a subseq.uent paper (reference 5) called 
such lines limiting lines. The velocity at the limiting line 
is never subsonic. However, the true characteristics of such 
limiting lines and their significance were not investigated 
by Tollmiou at that time. Recently Ringlob (reference 6) ob- 
tained another particular solution of isentropic irrotational 
flow in which the maximum velocity reached is approximately 
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twice the looel sound volocity, ffor this flow also a limit- 
ing lino appeared beyond which the flow cannot continuo, 
ffurthermorej he found the singular character of the limiting 
line, that is, infinite acceleration and infinite pressure 
gradient. Ton Zarman (reference 7, particularly pp. 351—356) 
demonstrated this fact for tho general two-dimensional flow^ 
He also suggested that the limiting line is the envelope of 
the Mach waves (fig, 1) and thus can occur only in a super- 
sonic region^ 3e also took its appearanc'e as "the oritorifln 
for hreakdown of isentropio irrotational flow, This general 
two— dimensional theory was established later by both Bingleb 
(reference 8) and lEollmien (reference 9)* Xollmien corrected 
some mistakes in Einglob's paper and, in addition, showed 
that the flow definitely cannot continue beyond the limiting 
line, aJho later fact introduced a "forbidden region" in the 
flow bounded by the limiting line, Ihis physical ab surd ity 
can be avoided only by relaxing the condition of irrotation— 
ality. But, as stated previously, for non— viscous fluids^ 
the transition from a flow without vorticity to that with 
vorticity can be accomplished only by shock vaves^ which at 
the same time also cause an increase in. the entropy^ 

However, before it can be concluded that the appearance 
of a limiting line, or the envelope of Mach waves, is the 
general condition for breakdown of is.entropic irrotational 
flow, it must be proved that the singular behavior of limit- 
ing lines is general and not limited to^ two— dimensional flow, 
Ihis is the purpose of the present paper, first the proporTy 
of liniting line in axially symmetric flow will be investigated 
in detail, Then the gonera.1 three— dimens ional problem will be 
sketched.^ !Ehese investigations confirm the results ~of Singlcb, 
Von barman, and Sollmien for these more general case's, 

5?herofore, by considering only the .gteady flow of non— 
Y.i^QQUs fluids, the criterion for breakdown of ieentropic 
irrotational flow is the appearance of a limiting line. Tor 
the actual motion of a solid body, however, the flow is neither 
steady nor non— viscous. Small disturbances always occur and 
almost all real fluids have appreciable viscosity. The small 
disturbances in the flow introduce the CLueation of stability. 
In other words,, the solution found for isentropic irrotational 
flow may be unstable even before the appearance of the limit- 
ing line, and tends to transform itself into a rotational flow 
involving shock waves at the slightest disturbance. If this 
is the caee, the criterion concerns not the limiting lino^ b'ut 
the stability limit, This pro'blem has yet to be solved. 
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The effect of viscosity will be limited to the ■boundary- 
layer if the pressure along the surface in the flow direction 
never increases too rapidly, tDhen outside the boundary layer 
the flow is Isentropic and irrotat ional. If the gradient of 
proBBuro is too large, the boundary layer will separate from 
tho surface. However, at low velocities such separation only 
widens the wake of the body and changes the pressure distri- 
bution over the body. But if the boundary layer separates 
at a point where the velocity outeid© the boundary layer is 
supersonic, additional effects may appear, The flow outside 
the boundary layer in this case can be regarded approximately 
as that of a solid body not of the original contour but of i\ 
now contour including tho "dead water" region created by tho 
separation. It is then immediately clear that the ideal 
isentropic irrotational flow around this new contour may havo 
a limiting line. Hence, the actual flow then must involve 
shock waves. In other words, the separation of the boundary 
layer in tho supers onid region may induce a shock wave and 
thus extend its influence far beyond the region of separation, 
Furthermore, the steep adverse pressure gradient across a 
shock wave may accentuate the separation, This interaction *- 
between the separation and the shook wave is frequently ob- 
served in experiments, 

She above considerations indicate the possibility of 
the breakdown of isentropic irrotational flow outside the 
boundary layer even before the appearance of the limiting 
line, Therefore, the Mach number of the undisturbed flow 
at which the limiting line appears may be called the "upper 
critical Mnch number," On the other hand, since shock waves 
can occur only in supersonic flow, the Mach number of the un- 
disturbed flow at which tho local velocity reaches tho veloc- 
ity of sound may be called the "lower critical Mach number,-" 
The actual critical Mach number for the appearance of shock 
waves and the compressibility burble must lie between these 
two limits. By carefully designing the contour of tho body ^ 
to avoid the crowding together of Mach waves to form an enve- 
lope and to eliminate adverse pressure gradients along the 
surface of the body, the compressibility burble can be delayed, 

AX I ALLY SYMMEa?EIG FLOW 



Iho solution of the exact differential equations for an 
axially symmetric isentropic irrotational flow was first given 
by Frankle (reference 10), The method was developed independ- 
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ently "by I'errari (reference 11), Their method applies par- 
ticularly to the case of Bupersonic flow ove? a "body of 
revolution, with pointed nose. In this case, the flow at 
the nose can "be approximated by the well-->known' solution 
for a cone, From this solution, the differential eq.uation 
is solved step "by step, using the net of characteristics 
which are real for supersonic velocities. In the following 
investigation, the chief concern is not the solution Of the 
partial differential oq.uation hut rather the occurrence and 
the properties of the limiting line' in an isentropic irrota— 
tional flov/. The general plan of attack is that of lollmien 
(reference 9), However, here the calculation is "based on 
the Legendre transformation of velocity potential instead of 
the stream function. 

If q. is the magnitude of the velocity, a the corre- 
sponding velocity of sound assuming an isentropic process, 
p the pressure, and p the density of fluid, the Bernoulli 
eq.uation gives 



a 1 — — — ~ •■' "A ' a I 1 + 



2 a, 



2 aa/ 



V 

8 \V-1 



2 ao^ / 



(2) 



(3) 



In those eq.uations, the subscript o denotes quantities 
corresponding to q = 0, and V Is the ratio of specific 
heats of the fluid. Let the axis of symmetry bo the x^axis, 
the distance normal to x— axis be denoted by y, and the 
velocity components^ along these two directions be denoted 
by u and v, respectively (fig, 2), The x-y plane Is, '■' 
theroforo, a meridian plane. Then the kinemat ical relations 
of the flow are given by the vorticlty equation 

- "tty = 0* (4) 



•Throughout this paper partial derivatives are denoted 
by subscripts. Thus "''x ^ t ^7 * "Ip 
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and the contiuuity equation . - v 

Eq.uatioiis (l) to (p), together with the relation .q.*au^+v^, 
specify the flow conipletely, 

Ho simplify the prohlem, a- velocity potential tp defined 
as follows is introducodt 



t 

« 

0?hen equation (4) is identically satisfied and equation (5), 
together with equations (l) and (2), gives the equation' for 9, 



She characteristics of this differential equation, to be 
called the characteristics in the physical plane, are given 
t>y g (x, y) = 0, where g (x, y) is determined hy the 
following equation 

It can "be easily seen from this equation that g is real 
only when q > a, therefore, the characteristics are real 
only in supersonic regions of the flow, 

The meaning of characteristics In the physical plane is 
immediately clear if one calculates the relation between the 
slope of a characteristic and the slope of a stream line in 
the meridian or x— y plane. By the definition of the function 
e(x,y), the value of g is zero, or constant,- along a char- 
acteristic, therefore, hy writing a quantity evaluated at a 
certain constant value of a parameter with that parameter as 
a subscript, the slope of the characteristic in the physical 
plane is 
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g 



gy 



(9) 



Along a streamline, the stream function defined by 

following cq.uations is constant; 

\lf =y-£-u, \lf s-,y-E-v 

^ Po ^ Po 

!Chorof oro, the slope of a streamline is 



(10) 



B^uatlons (8), (9), and (ll) give 




1 - ^ 
a^ 



± tan p 



(11) 



"•1 a 

vhere 3 is the Maoh angle given ty p= sin — . Therefore, 

equation (l2) shows that the characteristics in the physical 
plane are inclined to the streamlines by an angle equal to 
the Mach angle. Such linos are the wave froats of infinites- 
imal disturbances and aro called Mach waves. In other words, 
characteristics in physical planes are the Mach waves in 
that plane, There are two families of Mach waves inclined 
symmetrically with respect to each streamline. 

If to each pair of values of u and v, there is j2Jt§ 
pair of vPvlues of • x, y, then x and y can bo considered 
as functions of u, v. In other words, instead of taking 
X and y as independent variables u, v can be used as in- 
dependent variables, The plane with u and v as coordinates 
is called the "hodograph plane," An equation in tho hodograph 
piano corresponding to equation (?) can bo obtained by moans 
of Legendre*s transformation. By yrriting 
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X = + vy — 9 (is) 

it is soon thrt 

\ ' « y (14) 



Thon cq.uation (?) can "be written as 



(l-Z!)x,,fi[x„„X„-xj]-0 (16) 



She characteristics of equation ( 15 ) are given by f(u,v) s 0 
where f is the solution of the following differential eq.ua-* 

t ion 



Equation (16) shows that the characteristics in the hodograph 
plane depend upon the values of the derivatives of X which 
must Tao olDtained from equation (15), In other words, the 
characteristics in the hodograph plane change with the flow 
and arc not a constant set of curves as are those in two- 
dimensional prohlems. 

Ho obtain the relation between the characteristics in 
the physical plane and those in the hodograph plane, it is 
noticed that equation (9) can "be rewritten as 



(dy)g t (dx)g = - ; . (17) 
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n?hen eq,uatlOtt (8) is equivalent to 

However, in general, eq.uation (l4) gives the following 
relation between the differentials of x and y and those 
of u and v: 



(19) 

dy = X^^ du + dv 



into anog^uation for (^^)g ^^^"^s* '^^^ transformed 

equation can "be simplified by using ecLUation (l5), 2he final 



By means of these relations, equation (iS) can "be transformed 
into an og^ua'^ •* — f ^••'^ /j— \ «... 

equation ci 
relation is 



2 



- ' (f^i <--)g -{G-f^)-^ X^dv) ^ 



^ -V.J g 

(SO) 



= 0 



Therefore, if the first factor of equation (20) Is not zero 
the variations (<i^)g and (Av)^ along a characteristic in 

physical piano must' satisfy the relation 

{Q-T-y.t M ' - 1 ^ <^^>. 

0 (21) 
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Ihis is the same relation for the variations (du)f and 
(d.v)^ along a characteristic in the holograph plane as can 
"b© seen from eq.uation (l6) and the following relation ob- 
tained froEi the definition of f 

(dT)j : (du)j •= - * ^225 

!Dhe transformed characteristics of the physical plane and 
the characteristics of the hodograph plane themselves satisfy 
then the same first order differential equation, Iheref or " " 
these two types of curves are the same. In other words, the 
characteristics of the hpdograph plane are the representation 
of Maoh v^aves in the u—v plaiie» . 

TEE LIMITING: LISTE 
Equation (SO) shows that if 

X^u ^v *" ^v ~ ^ 



« - (22) 



then the transformed differential equation for the character- 
istics of the physical plane, or Kach vaves, is satisfied, 
;i?horeforo, if there is a line in the hodograph plane along 
which tho values of the derivatives of X are such that 
equation (23) is true, then this line when transferred to 
the physical plane will have its slope equal to that of one 
family of Kacsh waves. Such lines are called the limiting 
hodograph in u—v plane and the limiting line in physical 
plane. Since Mach waves occur only in the supersonic regions 
it is then evident that tho limiting line must appear in 
these regions. The significance of the adjective "limiting" 
will be made clear as other properties of such lines are 
investigated, 

How the question arises: Can the limiting hodograph 
be a characteristic in tt-v plane? Along a limiting hodo- 
graph, equation (23) gives 

^v ^ ^ _ ^uuu ^v ^ ^v ^uv ^ ^u \ivv ^ 2^ ^ 

\ Xuuv X.VV "* 2 Xtiv ^uvv + ^^u ^Wvv 
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where the subscript \ denotes the value along a limiting 
hodograph» ^ow the general differential eq.uation for 
equation (15), is true for the whole u— v plane* there- 
fore; the equation is still true hy differentiating it with 
respect to u and v. Ihe results can "oe simplified hy 
using equation (15) itself and equation (23), S?hen at the 
. limiting hodograph. 



+ JL Xuu Xuvv = (*'+!) XvT- 2 Xuv+ (V-1) -% >^ut (25a) 
Xv J a a* a 



+ 



rr ^ul ^vv» (Y-1) X^v- 2 C^ + l) ^u (25h) 

V J a® a* a* 



Equations (24), (25a), and (25h) are the only avgiiable 
equations involving no higher derivative than the third. On 
the othor hand, the slope of a characteristic in the hodo- 
graph plane can "be calculated hy equation (22), 

Shis equation together with equation (l6) gives 

Thercforo, if the limiting hodogrc^ph is a characteristic, then 

/J^ZN must satisfy equation (27), Hov.'evor, a simple calcula- 
Vdu/^ 
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tlon shows that it is not even possible to olitain a relation 
between (■^ ) and other q.uantities not involving the third— 

order derivatives of X, Hence, does not satisfy 

VduA 

equation (27). In other words, the limiting hodograph is 
not a characteristic. Transferred to the physical plane, 
this means that the limiting line is not a Maoh wave. But 
as shown in previous paragraphs, the limiting line is every- 
where tangent to one family of Mach waves, Conseq,uently, the 
limiting line must be the envelope of a family of Mach waves, 
fhis property of the limiting line can be taken as its physical 
definition. 



LXHia!IN& HODOGRAPH AND THE STESAMLINES 



At the limiting hodograph both eq.uations (15) and (23) 

hold. By eliminating one of the second— order derivatives, 
say X^y^i the following relation is obtained 

(Xvv)^ " CXuvh (28) 




a . 



The sign before the radical in eq.uation (28) can be either 
positive or negative, but not both, 5Phis relation will be 
used presently to show that the streamlines and one family 
of characteristics are tangent in the tt-v plane. 

From eq.uation (10), the differential of the stream 
function can be calculated as 



dvl/ = ~ y V dx + y JE- u dy (29) 

Po Po 

In this equation, y can be replaced by X^ according to 
equation (14) and the differentials dx and . dy replaced 
by the differential du and dv according to equation (19). 
Then — 
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Along a streamlino, 44/ = Oj therefore the slope of the stream- 
line in the hodograph plane is given by 

AA ^ V Xu^ - u Xuv ^^^^ 



At the limiting hodograph, equation (23) holds; therefore, 
equation (Sl) together with equation (28) gives 




u« 



v,,--UvA"^^,yarn 

where the sign before the radical oan be either negative or 
positive corresponding to the sign in equation (28), 

On the other hand, the slop© of the characteristics in 
the hodograph plane is determined by equation (27), By 

solving for (j^^ aad simplifying the result with the 
aid of equation (l5). 




Ihe sign before the radical is either positive or negative 
corresponding to the two families of characteristics. By 
using the positive sign in conjunction with the positive sign 
in equation (28), and similarly for the negative sign, 

a 
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Equations (32) and (34) show that the streamlines and on© 
family of characteristics are tangent to each, other at the 
limiting hodograph. This result is the same as that ob- 
tained for two— dimensional flow, (See references 7, 8, and 
9.) Those equations when compared with equation (l3) for 
the slope of Kach waves in the physical plahe yields the 
interesting result that the streamlines and one family of 
characteristics at the limiting hodograph are perpendicular 
to the corr-esponding Mach waves at the limiting line. 

Since 




(35) 



Equation (32) gives the following equation which holds at 
the limiting hodograph 

This equation can he reduced to more familiar form by intro- 
ducing the polar coordinates in the u-v plane: 



u = q cos e, V = q sin 9 



where 9 is the angle between the velocity vector and T>-axi 
Then equation (36) takes the form 

This can be regarded as the equivalent of equation (23) for 
defining the limiting hodograph. A similar relation exists 
for two-»dimensional flow, (See references 7, 8, and 9,) 

Along a streamline, the ratio between (dv)^ and 

(du)^ is given by equation (31), By substituting this 
ratio in equation (l9), the differential (dx)^ and 

(dy)^ along a streamline is given as 
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ufx X ^] 



— V X + "n X 



si 

\v J 



y (38) 



— V Y + 11 Y 

'niv ^VV 



At the limiting line, eq.ua.tion (23) is satisfied, !Ph.©n eq.ua- 
tion (38) shows that at the limiting line, the streamline has 
a singularity. Or, more plainly, Tdx)^ and (dy)^ at these 

points are infinitesimals of higher order than (du)y|, and 

(dv)^. By writing & for the distance measured along a 

streamline, equation (38) gives immediately 



V X + u X 
^ uv vv 



fx X - X 3] 

(_ uu vv uv J 



(39) 



Similarly, 



V X^^ u X^^ 



4 X X 
L UU vv 



uv J 



(40) 



therefore, at the limiting line, the acceleration along a 
streamline is infinitely large, Furthermore, since the 
pressure gradient (Pg)^ along a streamline is 



(PS)^ = -p«l<lg='-p[u (Ug)^ + V (Vg)^] 



(41) 



the pressure gradient at the limiting line is also Infinitely 
large. 
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Such infinite acceleration and pressure gradient lead 
one to suspect that the fluid is thrown hack at the limiting 
line. In othor words, the streamlines are douTjled 'back at 
this lino of singularity, lo investigate whether this is 
true, tho character of the relation - X^^^ = 0 

along a streamline has to "be determined. If the derivative 
of this oxpression along a streamline is not zero, then 

Xuu ^vv " ^v^ ^^-^ only a siaple zero at tho intersection 
of the limiting line and the streamline, Gonseq.uently, the 
differentials (dz)^ and (dy)^ will change sign hy passing 

through the limiting hodograph in tt-v plane along a stream- 
line. Hence, the streamlines will douhle back and form a 
cusp at the limiting line, Ihe derivative of x^u ^vv ""^uv^ 

along the streamline can he calculated with the aid of 3q,U£v- 
tion (30) 



rA(x X -X^)1=X X -2)^ X +X X 

tdu "u-ii vv uv ' uuu vv uv uuv uu uv 

"l^^uv ^v"* ^ \iv \ivv* \iu '^vvv^ (42) 



^ ^UU ^IV 



— V X +u .X 

uv vv 



The e3:prossion on the right of equation (42) cannot be reduced 
to zero by the available relrttions, which consist of equation 
(23), og.uation (15)., and differentiated forms of equation (15), 
S^hcrcfore, the expression concerned generally has only a simple 
zero at tho limiting hodograph and the stroamlines are doubled 
back at the limiting line. It will be shown later that there 
is no solution possible beyond the limiting line. Hence, the 
name limiting line. 



SSTTELOPE OP OHARAOTBEISa?IOS Iff HODOaSAPH PLAITE AND 
LIBTES OP OOffSIAira! VELOOiajy IN PHYSICAL PLANE 



Sinco the limiting line is the envelope of the Mach v/ares 
in the physical plane, i^t is interesting to see whether there 
is also an envelope for tho characteristics in the hodograph 
plane, The characteristics in the Ur-v plane are determined 
by equation (^6), The envelope to them can be found by elim- 



inating 




between equation (26) and the following equation 
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which is obtained, by equating to zero the partial derivative 

^dv 



of eq_uation. (26) with respect to result can "be 

simplified "by equation ^15), and then it hecomes simply 

1 — + ^ « (44) 

a* a* 

fhis is satisfied by either 

a = 0 (45) 

or 

u+vasa (46) 

Tho first condition, equation (45), when substituted Into 
equation (26) gives 

which shows that the circle of maximum velocity correspond- 
ing to a = 0, is the envelope of the characteristics in hodo:- 
graph piano, !l?he second condition, equation (46), is the 
spurious solution, sinoo generally the characteristic at 
q = a is not a tangent to the circle q = a. Hence a <= 0 
is the only envelope, r 

Tho linos of constant velocity in the hodograph piano 
are simply circles, Therefore 



'dv \ u 

V 



(48) 



By means of this relation and equation (19), the slope of the 
lines of constant velocity is given as 
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This eq.uation together with ectuation (30) gives the foirpwing 
interesting relation 




(5.0) 



In other words, a line of constant velocity in the physical 
plane is perpendicular to the streamline in the holograph 
plane at the corresponding points. 



lEBS LOST SOltJTIOJT 



Throughout the previous calculation, the possibility of 
using the Logendro transformation, is assumed, Shis req.uires 
that for each pair of values of u, v there is one and 
only one pair of values of y. However, it is not al— 

wayc tyuo. It is possihle to have a numher of points in 
the physical plane having the same value of u and v. If 
this is the case, then evidently it is impossible to .solve 
for X and y from the pair of functions u* u (s:,y), 
V s: V (x,y)« Mathematically, the situation is expressed by 
saying that the Jacobian Mu, v)/5(x,y) vanishes in the 
physical plane. Or 



^y - V = ° ^^^^ 

However, this is also the condition for a functional relation 
between u and v; for example, v can be expressed as a 
function of u. In other words, u and v are not independ- 
ent. Hence if a solution is "lost" or not included in the 
family of solutions allowing Legendre transformation, then for 
that solution, 

V = v(u) (52) 

It is seen that eq.uation (51) is then identically statisfied. 

By eliminating p from the continuity eq^uation, there 
is obtained 
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fhis equation can "be rewritten in the following form lay using 
equation (52) 

fhe vorticity equation, equation (4), can be expressed as 

- «y - 0 (55) 

from equations (54) and (55), it is possible to solve for 
Uj^ and Uy, The result is 

1(1-^^-)^ 2^^^ fl-I!^ C^yluy = -^il (56b) 
LV aS/ a3 du \ a^/. W/ J ^ y du 

By differentiating the first of equation (56) with respect 
to y, the second with rospoct to x, the following relar- 
tion is obtained by subtraction; 

s 

u + 1 « 0 (57) 



, S X 

du 



therefore 



or 



dv _ f(y) — X 
du y 



f(y) - X 
y = 

dv 
du 



(58) 



where f(y) is an undetermined function of y. However, 
equation (55) shows that for lines of constant values of u 
where du = u^^ i^^)^ + tiy (dy)^ = 0, 
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SO 



m — ■ . . ? •>■. . a constant (59) 



Eencn, linos of constant values of u and t are straight 
lines, fhis restriction reduces the function f(y) in 
oauation (58) to a numerical constant. Put f (y) =» E, eoLua— 
tion (58) is then 




(50) 



Ehcrefore lines of constant value of u and v are radial 
lines passing through the point x » K, Thus the lost solu- 
tion is nothing but the well— known solution for the flow 
over a conical surface, 

From eci.uation (59), it is seen, that lines of constant 
velooitjr are perpondicUiiar to the tangent of the u— v curve 
at the corresponding points. By substituting the value of 
1 

y from equation (57) into eq.viation (56a), a relation between 
u and V is ohtaine^^} 




(61) 



Ihis is the differential equation for determining the hodograph 
representing the flow over a cone. Figure 8 shows the hodograph 
for a cone of 30° semivortex angle and with a velocity at the 
surface of cone equal to 0,35 c, The maximum velocity is c — 
that is, the value of q corresponding to a » 0. Pigure 3 
is drawn from data given by Taylor and Maccoll (reference 12), 



It raay well be mentioned here that the lost solution for 
the axially symmetric flow is not limitsd to supersonic veloc- 
ity as is the jease for two-dimensional flow. In fact, Taylor 
and Kaccoll show that for small forward velocity of the cone, 
supersonic velocities occur only Just after the head shock 
wave. The velocity decreases as the surface of the cone is 
approached, finally, it becomes subsonic for points near the 
surface of the cone, Figure 4 shows a few examples taken from 
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their calculations (reference 12), The dotted curves in the 
figure are the Mach waves, The dotted straight lines are 
the 'boundaries "between the supersonic and the subsonic regions, 
furthermore, spark photographs of a conical shell in actual 
flight taken "by Maccoll (reference 13) do not indicate the 
presence of -shock waves in regions of flow vhere such transi- 
tion from supersonic to suhsonic velocities is expected. 
Therefore, at least for this particular type of flow, a 
smooth traJisition through sonic velocity actually takes 
placo. 



00NTI3JUATI0N Off SOLUTION BEYOND THE LIMITING LINE 



Since it is shown In a previous paragraph that the 
stroamlines are generally turned back at the limiting line, 
the CLUGstion arises* la it possible to continue the solution 
beyond the limiting line? Of course, there are two ways of 
continuing the solution: The new solution is Joined either 
smoothly to the given solution at the limiting line or with 
a discontinuity. As shown before, the limiting line Is the 
envolopo of one family of the Mach waves. Then at every 
point of this line its direction differs from that of the 
stroamlino by an angle equal to the Mach angle. But the 
Mach angle is not zero except at points where the velocity 
of fluid has reached the maximum velocity and the ratio 

^ ss 0, Therefore, the limiting line generally docs not 

q. 

coincide with the streamline, and the discontinuity at the 
junction of the solution at the limiting line cannot be that 
of a vortex sheet. The only other type of discontinuity is 
the shock wave. However, the angle between the limiting line 
and the flow direction is .eq.ual to the Mach angle. Then ac- 
cording to the result of the theory of shook waves, the dis- 
continuity across such a line vanishes. In other words, 
there cannot be a discontinuity at the limiting line. There- 
fore, it is impossible to Join a new solution at the limit- 
ing line with a discontinuity. 

As to the second possibility of Joining a new solution 
smoothly at the limiting line, it is seen that the flow beyond 
the limiting line must be irrotational and isentropic since 
the limiting line cannot be a shock wave. There are only two 
typos of isentropic irrotational flow? namely, one that allows 
the Legcndre transformation, and ono that does not, the "lost 
solution," Investigate the second alternative first. If the 
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solution "beyoni. the limiting line "belongs to the so— called 
lost solution, then since the junction at the limiting line 
must be smooth, the values of u and v at the limiting 
line must also satisfy equation (61), But the slope 

limiting line is given "by eq.uation (34). The 

secoiid derivative (^^"s^l ^^^^ then involve the fourth 

order derivatives of X, Besides these expressions, the 
available relations are equations (15), (S3), (25a), (SSh), 
and three more equations obtained "by differentiating equations 
(25) \ifith respoot to u and v. However, it is still impos- 
sible for *° satisfy an equation like equation (61) 

where no derivative of x appears. Hence, the limiting 
hodograph does not satisfy the equation for the lost solu- 
tion. In other words, the lost solution cannot be used to 
continue the flow beyond the limiting line, 

7he only remaining possibility is to continue the flow 
smoothly by another solution obtainable by a Eegendre trans- 
formation. Smooth continuation means that the values of u, 
V, and p must be the same at the Junction, the limiting 
line. Since shock waves do not appear, isentroplc relations 
still hold, Ihe density p is determined by velocity only, 
Ihe value of u and v are determined by the coordinates 
in the hodograph plane. The position of the limiting line 
In the physical plane is determined by X X Therefore, 

the problem can be stated as follows: At a certain given 
curve u(X), v(\) in the hodograph plane, the limiting 
hodograph, the values of X^, X^ are given; X is the 

parameter along the gi-ven curve. It is required to deter- 
mine a new solution of the differential equation (equation 
(15)) with these initial values, First of all, it is seen 
that with the given data, the left-«hand sides of the follow- 
ing equations are given: 



d\ 



(62a) 
(62b) 
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Therefore 




(63a) 
(63b) 



By substituting those values into equation (15), the second- 
degree terms reduce to 

/ 

\r " \v LdX dX ▼ d\ dX u J / \g^^J 

v/hioh is linear in Xv,-* -Sfherefore Xy^ can "be uniq.uely 

determined by eq.uation (15), In other words, with the given 
data, the second order derivatives of X at the given curve 
u(\), v(X) can "be determined uniq.uely, in spite of the fact 
that the differential eq.uation (l5) is of second degree, 
Friodrichs and Lewy (reference 13) have shown that under 
those ciroumstances , the function X within a region H 
(fig, 5) bounded by two characteristics and the given curve 
is uniquely determined except for an additional constant. 
Consequently there can be only one solution corresponding to 
the given data at the limiting hodograph. However, this solu- 
tion is the very ono which gives the reverse flow at limiting 
line. Therefore, it is impossible to continue the solution 
beyond the limiting line ovea by a Legendre transformation. 

Since all three alternatives fail to offer a way of con- 
tinuing the solution, the limiting line is truly an impossible 
boundary to cross. In other words, tho region beyond the 
limiting line is a forbidden region. This physical absurdity 
can be resolved only by the breakdown of isentropic irrota-" 
tional flow. 
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&ENBBAI> fHESiV-DIHZHSZONAIi FLOW 



Ehe methods used in previous sections for investigating 
the axlally sjmmetrio flow can "be easily extended to the 
general three-dimensional case. In the present section, this 
investigation will "be sketched lirlefly and the results indi<- 
cated* 

Let the three components of velocity along the three 
coordinate axes x^ y, and z he denoted by u, v, and 
V, respectively, fhen ly introducing a velocity potential 
cp defined by' 

u = cp,., ^ " <Py» w " (65) 

the differential eg.uation for cp of an isentropic irrotational 
flow can he written as (reference 7) 



^^xx ^yy T^zz' 



^ ^xx"*"^ q>yy+ w <p^^ + 2vw<pyjj+ 2wu923, + 2uvcp3,y (66) 

If, for every triad of u, v, w, there is only one triad 
of X, y, z, then the Legendre transformation can he 
used. Xhus 

X = 11X + vy + wz — «p (67) 

and 

■^u ° ^» '^■v " y» ^v" ' (68) 

The differential eq.uation for qj, equation (66), is then 

transformed into 



a^. bO-ff'^+OA-G +AB-h''1= ^'^ (BO-]P^) + v^ (OA-G^) 

8 8 ^^^^ 

+ w <AB- S ) + 2vw (GE-^ AI) + 2wu (HP- BG) + 2uv (IQ-OH) 
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where the following notations are used 



0 = X 



WW 



, J^X , --^ , 



vw' wu" 



H = X 



(70) 



By analogy with the axially symmetric case, the limit- 
ing hodograph surface is defined as the surface in the u, 
V, w space, or hodo'graph space, where the following rela- 
tion holdsS 



A E a 

H B I" 
G ff 0 



(71) 



The properties of this limiting hodograph and the correspond- 
ing limiting surface can be found ty considering the behaTior 
of stroamlineB and characteristics at such surfaces. 

I'rora eg^uation (68) the differentials of x, y, and z 
can he written as 



dx a A du + H dv + & dw 



dy B H du + B dv + I" dw 



da «s G> du + r dv + 0 dw 



(73a) 
(72h) 
(72c) 



Along a streamline, the differentials dx, dy, and dz must 
he proportional to u, v, and w, respectively, Ihus the 
eq.uation of a streamline in physical space is 



u 



(73) 



where the suhsoript ^ indicates values taken along the 
streamline, The eq.uation of a streamline in hodograph space 
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is obtained by eliminating dx, dy and dz from equation 
(73) "by equation (72). The result is 

^^^V <^^)m; 



(74) 



au + hv + gw hu + + f w gu + f v + cw 



where a Is the_co— factor of A in the' determinant A of 
equation (71), b the co— factor of B, and so forth. Equa- 
tion (74) can be used, in turn, to eliminate two of the three 
differentials du, dv, and dw in the right of equation 
(72)« The result is 

/ , V UA du ' / nc \ 

C4x)^ = ::: f'^a) 

^ au + hv + g w 

vA dv 

(<iy) = = r: =- (75b) 

* hu + bv + f w 



. . wA dw 

(dz) = (75o) 

* gu + fv + cw 

At the limiting surfaoe, A = 0 as defined by equation (71)' 
therefore the streamlines have a singularity there. Similar 
to the axially symmetric flow, the streamlines generally ard 
turned back and form a cusp at this surface. The accelera- 
tion and the pressure gradient are, of course, infinitely 
large at such places. 

The characteristic surface g(x,y,2) =0 in physical 
space is determined by the equation 



+ gy + gz J = 



^Sg^s^. T2ey3+ v'^s^^+ 2vwgygj^+ 2wug^g^+ Suvg^gy (76) 
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Since this eq.uatioii is a second— degree equation, there are 
two families of surfaces passing through each point, Eheso 
surfaces are the wave fronts of infinitesimal disturTjanoea 
in the flow and can "be called the Mach surfaces, The char—, 
actoristic surface f(u, v, w) = 0 in the hodograph space 
is determined by the equation 

j^(^0)f^^+ (0+A)f/+ (A+B)f^^- Sff f^f^- 2& f^f^- 2H f^'f^J 
- 3" f A f„ f„1 + 2wu fsff+Pff - &f^-Bffl 

U VW VW UV V wu 

+ 2UV ^Qr f^f^ + y f^f^ - H f^= - 0 f^ j (77) 

By transforming equation (76) for Maoh surfaces to hodograph 
space, it can he shown that the tran&formed equation is 
satisfied either "by tlie characteristics in hodograph space 
determined "by equation (77) or hy the limiting hodograph 
determined "by equation (71)« Therefore, here again the 
limiting' surface is the envelope of a family of Mach surfaces. 

By using equations (74) and (77), it is possible to show 
that the streamlines In the hodograph space are tangent to 
the characteristic surfaces at the limiting hodograph, Fur- 
thermore, by using equations (69), (71), and (74), the in- 
clination of the streamlines at the limiting hodograph can 
be calculated. In fact, if (ds)^ = ( du) ® + ( dv) ^ + ( dw) ^, 

q® = u'^ + v'^ + w'^, the following relation is obtained 



= ^ or - - (78) 



This relation Is really equivalent to equation (32), In 
other words, at the limiting hodograph, the inclination of 
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the streamlines and characteristics to the q. = constant 
Burfaco is eq.ual to the Mach angle (fig, 6), It thus seems 
the "breakdown of general steady isentropic irrotational flow 
of nonviscous fluid is connected with the appearance of the 
envelope of Mach waves in physical space and the tangency of 
streamlines and characteristics in hodograph space, 

California Institute of Technology, 

Pasadena, Calif., August 24, 1943. 
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figure 1.- Limiting line as the 

envelope of Uach waves. 
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Figure 2.- Stream line and velocity 
components in an axially 
symmetric flow. 
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PHrS/CAL PLANE 

Figure 3.- Hodograph of the flow over a cone of 30° half vertex angle 
and a surface velocity q enual to 0.350. 
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Figure 4.- Flow oyer cones of various vertex angles involving 
rr^^r. subsonic rsgions . eg = half vertex aHgle, us » 

velocity over the surface of cone. • s 




PigurB 5.- Region R whsrs tte solution Is uniquely 
deteroined by given data at the Uniting 

bodograpb. 
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Figuxe 6.- Gdometrl^al relations between streaa line 

4nd cbariioterlatlc surface at limiting 
Burfaoe la nodograp:i space. 
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